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We compute numerically the quasinormal modes of Kerr-Newman black holes in the scalar case, 
for which the perturbation equations are separable. Then we study different approximations to 
decouple electromagnetic and gravitational perturbations of the Kerr-Newman metric, computing 
the corresponding quasinormal modes. Our results suggest that the Teukolsky-like equation derived 
by Dudley and Finley gives a good approximation to the dynamics of a rotating charged black hole 
for Q < M/2. Though insufficient to deal with Kerr-Newman based models of elementary particles, 
the Dudley-Finley equation should be adequate for astrophysical applications. 
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I. INTRODUCTION 



The electrovacuum black hole solutions of the Einstein-Maxwell system can be uniquely described by the Kerr- 
\ Newman (KN) metric pj, which is the most general among classical black hole solutions. The KN metric is specified 
■ by three parameters: the black hole mass M, the charge Q and the angular momentum per unit mass a = J/M. As 
long as M 2 > Q 2 + a 2 the KN metric describes a black hole, otherwise it has a naked ring-like singularity. As Q — > 
the KN metric reduces to the rotating Kerr metric, and as a — > it reduces to the Reissner-Nordstrom (RN) metric. 
Both limits have been studied in great detail 0, • 

As early as 1968 Carter realized that the KN solution has a magnetic dipole moment corresponding to the same 
g-factor g — 2 as the electron Q. This led to the suggestion (recently revisited in |f|: see discussion and references 
therein) that the KN metric could provide a reasonably adequate model for the external Einstein-Maxwell field of 
elementary particles 0, Q . 

For astrophysical black holes Q is likely to be negligible, electric charge being shorted out by the surrounding 
I , plasma @- If any charge is present, the no-hair theorem guarantees that charged rotating astrophysical black holes 
^Jy are described by the KN metric. Punsly proposed a gamma-ray burst model based on bipolar outflow from a fast- 
rotating black hole endowed with a small (by gravitational standards) charge and surrounded by a magnetosphere . 
Van Putten suggested that hypernovae or black hole-neutron star coalescence may generate rapidly spinning black 
hole-torus systems, and that the spin energy of the hole could power gamma-ray bursts [l^. In his model the black 
hole-torus system has a gravitationally weak magnetic field as a result of the remnant flux from the progenitor star 
(a massive star in hypernovae or a neutron star for a coalescence remnant), and the interplay between rotational 
effects and magnetic fields may be relevant. Ruffini et al. noticed that pair creation induced by the Heisenberg-Euler- 
Schwinger vacuum polarization during collapse leading to formation of a charged black hole could explain the energetics 
of gamma-ray bursts, and proposed a detailed collapse model in a series of papers More recently, Araya-Gochez 
suggested that intermittent hyper-accretion produced by magnetorotational instabilities in the accretion disc of a 
rapidly rotating, newly born black hole may induce resonant excitation of the black hole's quasinormal modes. His 
estimates suggest that a 15 black hole spinning at a ~ 0.98M and located at 27 Mpc could produce gravitational 
waves detectable by LIGO II [l3 ■ 

Both in elementary particle models and in astrophysical applications rotation plays a crucial role. In geometrical 
units, the typical angular momentum of an electron a ~ K/2m e — 1.93 x 1CP 11 cm, the length associated to the electron 
charge Q — e — 1.38 x 10~ 34 cm and the corresponding mass scale M — m e = 6.76 x 10~ 56 cm, so that M cQCa 
and the gravitational field is spin-dominated. In fact, the rotation parameter is so large that the KN metric can 
only model the Einstein-Maxwell field of the electron outside some small radius ro surrounding the ring-like naked 
singularity. Astrophysical black holes are also expected to be formed in rapid rotation: recent simulations suggest 
that supermassive stars would form black holes with a/M ~ 0.75, and many observations are consistent with black 
holes spinning close to the extremal limit a/M ~ 1. Reference [13J provides an updated discussion of the observational 
evidence for black hole spins and of the related uncertainties in present-day astrophysical models. Most black hole 
models for gamma-ray bursts require rapid rotation 0, llOj , and the magnetorotational instability considered in [l2| 
is most efficient for large spin parameters, a/M > 0.9. 
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The KN solution is the only asymptotically flat solution of the Einstein-Maxwell system for which the geodesic 
and Klein-Gordon equations can be solved by separation of variables |14j|. The Dirac equation in the KN metric 
is also known to be separable ^]|. Scalar and Dirac perturbations of a KN black hole can therefore be treated 
using the same general methods that apply to Kerr black holes. In particular, it is straightforward to compute the 
quasinormal modes (QNMs) of scalar perturbations of the KN black hole. As far as we know, a complete analysis 
of the corresponding QNM spectrum is still lacking. The first objective of this paper is to fill this gap, presenting a 
complete continued-fraction calculation of scalar QNMs of the KN metric for all values of a and Q. 

Studies of the interplay of electromagnetic (EM) and gravitational fields in the KN metric are plagued by a major 
technical difficulty: all attempts to decouple the EM and gravitational perturbations of the KN spacetime to date 
have failed. Section 111 of gives an introduction to this long-standing unsolved problem. Dudley and Finley (0|, 
henceforth DF; see also 01) make a remarkable study of the separability of linear perturbations of the solutions of the 
Einstein-Maxwell equations found by Plebahski and Demiahski . I'll , which include all vacuum Type D solutions . 
In their work Dudley and Finley "either keep the geometry fixed and perturb the electric field or, of more interest, 
keep the electric field fixed and perturb the geometry" [l|| . This approach should be appropriate for values of the 
charge Q at most as large as the perturbations of the spacetime metric (in geometrized units). Dudley and Finley 
show that a sufficient condition for decoupling is that the spacetime be of Type D, and that the decoupled equations 
only separate (in Plebahski-Demiahski coordinates) for perturbing fields of spin s — 0, 1/2, 1 and 2. 



Mashhoon |20| first presented arguments in favor of the stability of the KN metric. Instead of explicitly computing 
QNMs, he used an approximate argument (originally due to Goebel [21| and reviewed in the Appendix): in this sense 
his stability proof is not fully convincing. Mashhoon's analysis is based on perturbations of test null rays in the 
unstable circular orbit of a KN black hole, and is strictly valid only in the eikonal limit / ^> 1. Here we show by a 
direct calculation that Mashhoon's claim is correct, and that his predictions are surprisingly accurate even for small 
values of I. At least for scalar perturbations our investigation of the QNM spectrum can be considered conclusive, 
because for s = the DF equation is exact. 

One of us (KK) used the DF equation to compute the fundamental gravitational QNM using WKB methods [22] ■ 
The main problem of this approach is not computational, but physical. The WKB approximation is reasonably 
accurate for all values of a and Q (as we show comparing results to a continued- fraction calculation in Table ITTf> . 
However, the DF equation is derived under a number of mathematical assumptions, and it is only approximately valid 
for gravitational and EM perturbations of KN black holes. 

A purpose of this paper is to clarify the physical range of validity of the DF equation and the physical meaning 
of their approximations. We first notice that the a — > limit of the DF equation does not yield any of the two 
Schrodinger-like equations describing coupled EM-gravitational perturbations of the RN black hole. Then we consider 
metric perturbations of the RN black hole freezing EM perturbations (or, vice versa, EM perturbations freezing the 
metric), and we compute the associated QNMs. We show that for values of the charge Q < M/2 the results are 
in good (but not exact) agreement with the DF equation. In other words, the assumptions behind the separability 
conditions leading to the DF equation are not equivalent to simply freezing EM (or gravitational) perturbations. 
Nonetheless, the DF equation yields EM and gravitational QNMs in good quantitative agreement with results for 
the coupled EM-gravitational RN perturbation equations when Q < M/2. Furthermore, the qualitative behavior 
of gravitational and EM QNMs as functions of charge and angular momentum is very similar to the results of our 
"exact" calculations for scalar perturbations. Our calculations suggest that the DF equation provides a reasonable 
approximation of the EM and gravitational dynamics of a KN black hole, at least for the gravitationally small values 
of Q expected in astrophysical applications. 

The plan of the paper is as follows. In Sec. ^ we present the DF equation and we compute the corresponding 
QNM frequencies using Leaver's continued fraction technique [23 • Our calculation is "exact" (in the sense that the 
perturbation equations involve no approximations) for s = 0: as far as we know, it represents the first example of 
numerical evidence for the stability of the KN metric to scalar perturbations. For EM and gravitational perturbations 
the RN (a — > 0) limit of the DF equation does not yield the standard RN perturbation equations, which are summarized 
for completeness in Sec. 11111 To clarify the meaning (and the limits of applicability) of the DF equation for EM and 
gravitational perturbations, in Sec. IIVI we consider approximate versions of the RN perturbation equations. To our 
knowledge, the "frozen" picture we present in this Section has never been studied before. To conclude we summarize 
our findings, discuss their meaning for concrete applications in astrophysics and elementary particle models, and list 
some open problems for future research. In the Appendix, for completeness, we summarize results of Ref. |20j on 
perturbed equatorial circular orbits of null rays in the KN metric and their relation with the QNMs in the eikonal 
limit. 
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II. DUDLEY-FINLEY PERTURBATION EQUATIONS FOR KERR-NEWMAN BLACK HOLES 



Dudley and Finlcy reduced the KN perturbation problem to a pair of differential equations: one for the angular 
part of the perturbations, and the other for the radial part 0, In this paper we consistently adopt Leaver's 
conventions |2^|. In particular, we use geometrical units and set 2M = 1. In Boyer-Lindquist coordinates, defining 
u — cos #, the angular equation 
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is the same as in the Kerr case. The angular separation constant E used in 22] is related to Leaver's Ai m by 
E = Ai m + s(s + 1). Correcting a typo in Eq. (5) of (the factor A\sK' should read 2isK'), the DF radial equation 
reads 
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where K = (r 2 + a 2 )u — am, A = r 2 — r + a 2 + Q 2 and A = Ai m + (aui) 2 — 2amui. The parameter s — 0, —1, —2 
for scalar, EM and gravitational perturbations respectively, and a is the Kerr rotation parameter (0 < a < 1/2). In 
the Schwarzschild limit a — > the angular separation constant A[ m — > 1(1 + 1) — s(s + 1). The DF equation was also 
studied by Detweiler and Ove j^. It is exact only when the spin s = or when we set the charge Q — (and then 
it reduces to the standard Teukolsky equation for Kerr black holes |25j). 

Let us consider the boundary conditions of Eq. J2J at the horizon. The horizon radius r+ = (1 + b)/2, where 
b = yl — 4(a 2 + Q 2 ). From the indicial equation it follows that the ingoing solution at the horizon is such that 
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where <r+ = [w(^+ — Q 2 ) — am]/b. Similarly, imposing purely outgoing radiation at infinity we find 



(4) 



Both these boundary conditions on the radial equation J5J and the corresponding regularity conditions on the angular 
equation can be cast as three-term continued fraction relations of the form 
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The coefficients of the radial continued fraction are: 



= f + (2 - 2ia+ - s)j (6) 
+ {— 32(wa 3 m + i6 2 cr + ) + 16[i& 2 ser + + isamb + uj 2 a 4 + a 2 m 2 + a 2 uj 2 b + b 2 — b 2 a 2 _ — amtub — \suja 2 b — sb 2 ] 
+ 8[b 2 a 2 uj 2 + a 2 uj 2 - amuj - ib 2 suj - b 2 amuj] + A[cu 2 b + uj 2 b 3 ~ ib 3 suj - isub] + u 2 (b 4 + 66 2 + 1)}/(166 2 ) , 



Pi = -2j 2 + 2Mfo+l) + 2iff + -l]j 

+ {32wa 3 m + 16[i6 2 cr+ + ujb 3 a + + b 2 a+uj + b 2 a\ 



+ 8{b 2 amu) + iujb 3 +uj 2 b 3 + ib 2 uj + amu; - 
= j 2 - [2i(u + a +)- s]j 



sb 2 



2 2 

a oj 



2 2 2 4l 

am — u a \ 

b 2 -b 2 A lm ]+u 2 {3b 4 



6b 2 



l)}/(86 2 



{32{oja 3 m + b 2 a + uj] + 24i6 2 scij + 16[6 2 cr^ + ib 2 s<j + + a 2 ui 2 b + isamb — amojb — isuua 2 b — uj 2 a A 



2 2l 

a m 



+ 10w 2 6 2 + 8[6 2 



b 2 a 2 w 2 - a 2 u} 2 } + A[uj 2 b 3 + uj 2 b - ib 3 su - isub] - uj 2 (1 + 6 4 )}/(166 2 ) . 



and those of the angular continued fraction can be found in [23|. To find QNM frequencies we first fix the values of 
a, £, m and w, and find the angular separation constant A; m (w) looking for zeros of the angular continued fraction. 
Then we use the corresponding eigenvalue to look for zeros of the radial continued fraction as a function of lo. The 
n— th quasinormal frequency is (numerically) the most stable root of the n— th inversion of the continued-fraction 
relation JHJ, i.e., it is the root of 
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TABLE I: Fundamental scalar QNM of the KN metric with I — m — for selected values of the charge and angular momentum. 



_o Q = 

(0.220910,-0.209791) 

0.1 (0.221535,-0.209025) 

0.2 (0.223398,-0.206506) 

0.3 (0.226342,-0.201397) 

0.4 (0.229074,-0.191402) 



Q = 0.1 
(0.222479,-0.210122) 
(0.223121,-0.209319) 
(0.225033,-0.206673) 
(0.228021,-0.201267) 
(0.230486,-0.190541) 



Q = 0.2 

(0.227471,-0.210998) 
(0.228169,-0.210067) 
(0.230230,-0.206961) 
(0.233297,-0.200439) 
(0.234035,-0.186969) 



Q = 0.3 
(0.236910,-0.211849) 
(0.237712,-0.210596) 
(0.239994,-0.206302) 
(0.242568,-0.196687) 



Q = 0.4 
(0.253168,-0.210094) 
(0.254019,-0.207839) 
(0.255542,-0.199555) 



The infinite continued fraction appearing in equation JJJ can be summed "bottom to top" starting from some large 
truncation index N. Nollert |26| has shown that the convergence of the procedure improves if such a sum is started 
using a series expansion for the "rest" of the continued fraction, i?jv, defined by the equation 

Rn = - r 1N+1 ff ■ (8) 

The series expansion reads 

oo 
k=0 

where the first few coefficients have the same form as in the Kerr case [3(j , except for the charge-dependent correction 
in b: C = -1, C 1 = ±V^2kJ&, C 2 = [3/4 + iuj(b + 1) - a]. 

In Tabled we give numerical results for the fundamental scalar QNM with I = m = as a function of charge and 
angular momentum. Dashed entries in this and the following Tables correspond to combinations of Q and a for which 
Q 2 + a 2 > M 2 . 

Fig- El shows trajectories described by selected scalar QNMs in the complex plane. In each panel, the thick black 
line corresponds to modes of a RN (possibly charged but non-rotating) black hole. As Q increases the RN mode moves 
counterclockwise in the complex frequency plane; open circles on the thick line mark increasing values of the charge 
(Q = 0, 0.05, 0.1, . . .). For five fixed values of Q (namely Q = 0, 0.1, 0.2, 0.3, 0.4) we plot the KN QNM trajectories 
as we increase a: the results are the five thin lines branching from the RN limit. Open circles on these curves mark 
increasing values of the angular momentum (a = 0, 0.05, 0.1, . . .). The observed spiralingbehavior (as a function of 
both charge Q and angular momentum a) confirms and extends results presented in [23, [23, l29l l3fl l3lt 1321 l33l | and 
summarized in |34| . The bottom-right panel shows trajectories of the fundamental scalar mode with I = 2, m — 0. 
As we increase I the counterclockwise bending (that for I = is very pronounced even for the lowest overtones) only 
shows up at higher overtone indices (n ~ 10 for 1 = 2). 

When I > 0, and for any value of the charge < Q < 1/2, rotation induces a Zeeman-like splitting of the modes. 
In Fig. |21 we show the splitting of the real part of the frequency of the fundamental scalar QNM with I = 2 for 
three representative values of the charge (from left to right: Q — 0, 0.2, 0.4). For Kerr black holes, the lowest-lying 
modes with I = m tend to approach the critical frequency for superradiance mf2 (where f2 is the rotational velocity 
of the black hole horizon) in the extremal limit. Detweiler (3j| first presented analytical arguments to explain this 
clustering of modes at the superradiant frequency in the extremal limit. The implicit assumptions in his argument 
have recently been re-examined in |3(|[33. In t he g eneral KN case the extremal frequencies for modes with I = \m\ 
have analytically been computed by Mashhoon [20||. We quote his result, which is in excellent agreement with our 
numerical calculations (see Fig. [21 below), at the end of the Appendix. 

In Table ITTI we consider gravitational perturbations with I = 2, m = 0. To assess the reliability of WKB methods in 
the KN case we compare results from Leaver's continued fraction approach (which can be considered exact, within the 
given numerical accuracy) with the third order WKB technique used in [22|, building on previous results in |38, 39]. 
The agreement for the fundamental gravitational QNM with I = 2, m = is excellent, typical errors being smaller 
than one part in a thousand. The WKB approach systematically underestimates the real parts and overestimates 
the imaginary parts. Third order WKB methods become slightly less reliable for large rotation parameter a; at fixed 
a, their accuracy depends very weakly on Q. Table ITTTl shows results for the I = m = 2 "barlike" mode, which is 
gen erally believed to be dominant in the quasinormal ringing following gravitational collapse to a rotating black hole 

Perhaps the most important outcome of our analysis is a null result: we could not find any unstable mode. Our 
quasinormal mode calculation confirms qualitatively Mashhoon's arguments in favor of the stability of the KN metric 
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FIG. 1: Left to right and top to bottom, the first three panels show trajectories in the complex plane of the first three scalar 
QNMs of the KN black hole for I = m = 0. The thick black curve corresponds to a = (RN limit). The thin curves are 
obtained increasing a from zero to the extremal limit for fixed values of the charge (Q = 0, 0.1, 0.2, 0.3, 0.4 respectively). 
Open circles (when present) mark selected values of the charge (Q = 0, 0.05, 0.1 . . . along the thick RN line) and angular 
momentum (a = 0, 0.05, 0.1 .. . along the thin lines). Bottom right panel: same plot for the fundamental mode with s = 
and I — 2. Thin lines are now trajectories of the m = KN modes for Q — 0, 0.1, 0.2, 0.3, 0.4. For clarity we do not show 
modes with other values of m. 




FIG. 2: Real part of the fundamental scalar QNM of a KN black hole for I — 2 and different values of m. Curves from top to 
bottom refer to m = 2, 1, 0, — 1, —2. The three panels correspond to different values of Q, as indicated. For Q ^ 0, the vertical 
line marks the extremal limit. 
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FIG. 3: Real part (top) and imaginary part (bottom) of the fundamental I = 2 gravitational QNM of a KN black hole. Solid 
lines refer to m = 2, dashed lines to m = —2. Different panels correspond to different values of the charge, as indicated. For 
Q 7^ 0, the vertical line marks the extremal limit. Thick (black) lines are numerical results from the DF equation. Thin (red) 
lines are obtained from Mashhoon's arguments, which are strictly valid only as I — > oo (see Appendix). Considering that we are 
looking at modes with I = 2 the agreement of Mashhoon's predictions with numerical results is quite impressive, particularly 
for the real part of the fundamental QNM frequency and for large values of Q and a. 



TABLE II: Comparison of QNM frequencies of the DF equation obtained by the continued fraction method (first row for each 
value of a) with third-order WKB results from |22T | (second row). Numbers refer to the fundamental gravitational QNM with 
1 = 2, m = 0. 



a 


Q 


= 


Q = 


= 0.1 Q = 0.2 Q = 0.3 Q = 0.4 





(0.747343, 
(0.7463, 


-0.177925) 
-0.1784) 


(0.753643, 
(0.7530, 


-0.178602) (0.773902,-0.180620) (0.813278,-0.183793) (0.886468,-0.186603) 
-0.1792) (0.7734,-0.1812) (0.8128,-0.1844) (0.8860,-0.1870) 


0.1 


(0.750248, 
(0.7496, 


-0.177401) 
-0.1780) 


(0.756664, 
(0.7560, 


-0.178053) (0.777325,-0.179977) (0.817630,-0.182905) (0.893230,-0.184893) 
-0.1786) (0.7835,-0.1804) (0.8172,-0.1834) (0.8928,-0.1854) 


0.2 


(0.759363, 
(0.7586, 


-0.175653) 
-0.1762) 


(0.766159, 
(0.7654, 


-0.176211) (0.788148,-0.177790) (0.831579,-0.179789) (0.915681,-0.178292) 
-0.1766) (0.7874,-0.1782) (0.8310,-0.1802) (0.9150,-0.1788) 


0.275 (0.771072, 
(0.7702, 


-0.173158) 
-0.1736) 


(0.778394, 
(0.7776, 


-0.173563) (0.802246,-0.174561) (0.850222,-0.174875) (0.947708,-0.165369) 
-0.1740) (0.8014,-0.1750) (0.8494,-0.1754) (0.9470,-0.1660) 


0.3 


(0.776108, 
(0.7752, 


-0.171989) 
-0.1724) 


(0.783669, 
(0.7828, 


-0.172315) (0.808382,-0.173003) (0.858524,-0.172355) 
-0.1728) (0.8074,-0.1734) (0.8578,-0.1728) 


0.4 


(0.803835, 
(0.8026, 


-0.164313) 
-0.1648) 


(0.812886, 
(0.8118, 


-0.163977) (0.843068,-0.161921) 
-0.1646) (0.8410,-0.1626) 


0.45 


(0.824009, 
(0.8228, 


-0.156965) 
-0.1576) 


(0.834319, 
(0.8330, 


-0.155756) (0.869094,-0.149647) 
-0.1564) (0.8680,-0.1506) 


0.48 


(0.838981, 
(0.8378, 


-0.150055) 
-0.1510) 


(0.850244, 
(0.8490, 


-0.147806) 
-0.1488) 
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TABLE III: QNM frequencies of the DF equation obtained by the continued fraction method. Numbers refer to the fundamental 
gravitational QNM with I = m = 2. 



a Q = Q = 0.1 Q = 0.2 Q = 0.3 Q = 0.4 

(0.747343,-0.177925) (0.753643,-0.178602) (0.773902,-0.180620) (0.813278,-0.183793) (0.886468,-0.186603) 

0.1 (0.804291,-0.176622) (0.812131,-0.177253) (0.837615,-0.179080) (0.888479,-0.181653) (0.989234,-0.181863) 

0.2 (0.879684,-0.173764) (0.890062,-0.174238) (0.924394,-0.175451) (0.996178,-0.176132) (1.158394,-0.167066) 

0.275 (0.956470,-0.169604) (0.970151,-0.169774) (1.016478,-0.169756) (1.120103,-0.166278) (1.432934,-0.122011) 

0.3 (0.988090,-0.167530) (1.003384,-0.167519) (1.055787,-0.166749) (1.177413,-0.160442) 

0.4 (1.172034,-0.151259) (1.200883,-0.149228) (1.311304,-0.138796) 

0.45 (1.343229,-0.129738) (1.395907,-0.123109) (1.678651,-0.074535) 

0.48 (1.535348,-0.098867) (1.653330,-0.077666) 



|2Cj . Fig. |31 shows that the agreement is more than qualitative. In fact his argument (which is based on the eikonal 
approximation, and is expected to be accurate only for I 3> 1) captures almost perfectly the behavior of the funda- 
mental gravitational QNMs with I = \m\ = 2 as functions of charge Q and angular momentum a. This agreement is 
particularly good for the real part of the frequencies, and it can be exploited for "quick and dirty" estimates of QNM 
excitation in astrophysical scenarios (see e.g. [12 )■ 

A major problem of the DF equation is that it only provides an approximation to the problem of coupled EM- 
gravitational perturbations of the KN black hole. One can easily check that the limit Q — > of Eq. @ yields the 
Teukolsky equation [25J and the corresponding QNM frequencies. In fact, the first column of our Tablc[H]is in perfect 
agreement with Table 2 in |23|. However, as a — ► (first line of Table ITTf we do not recover the QNM frequencies 
of the RN metric. This is no surprise: for the reasons discussed in the Introduction we expect the DF equation to 
be a good approximation only for small values of Q. Below we will show that this is indeed the case, and give a 
quantitative meaning to the above statement. 



III. REISSNER-NORDSTROM: THE STANDARD TREATMENT 



Here we briefly summarize the computational procedure for RN black holes (a = 0). More details can be found in 
Refs. [23L l2(i| and especially [22 • This Section is primarily intended to establish notation for the "frozen" approxima- 
tion to be introduced in Section HVI 

It is well known that polar perturbations of the RN black hole can be obtained from the axial perturbations by 
a Chandrasekhar transformation [2; in particular, the QNM spectra of polar and axial perturbations are the same. 
For this reason in the following we only consider axial perturbations, that we denote by a superscript (— ). Let us 
introduce a tortoise coordinate by the usual relation 



dr _ A 
dr* r 2 



(10) 



where A = r 2 — r + Q 2 (recall that in our units 2M — 1 and < Q < 1/2). Explicitly, the tortoise coordinate can be 
written as 



^2 2 

— — — ln(r — r+) ln(r — r_) 

_ — r_ r + — r_ 



(11) 



where r± = (1 ± yl — 4Q 2 )/2 is the location of the inner (Cauchy) and outer (event) horizons of the RN metric. 
Using tensor spherical harmonics to separate the angular dependence and a Fourier decomposition to get rid of the 
time dependence, EM and axial gravitational perturbations of the RN metric are described by two coupled wave 
equations: 



drl 
dr 2 



H, 



(-) 



H 



(-) 



A 
A 



l{l + l)r 



1(1 



4Q 2 



4Q 2 



l)r-| 
' 2 r 



H. 



(-) 



\Ht ) +2Q y /(l- l){l + 2)H[ \ 



H 



(-) 



(12) 



+ -H[ ] +2Q^{l-l){l + 2)H, 



(-) 
2 
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where ' corresponds to perturbations of the gravitational field and H[ to perturbations of the EM field 0, ^2 • 
The usual procedure is to decouple the system (If 2fl to obtain Schrodinger-like equations of the form 



(13) 



where 



AO 2 

1(1 + l)r -qj + ^y- 



= 1,2, i^ j) , 



(14) 



and 



9i = 



3+ y/9 + l6Q 2 (l~l)(l + 2) /2 



'72 



3- y/9 + l6Q 2 (l-l)(l + 2) j2 



(15) 



The price to pay is that, unlike H[ ' and ' ' , the decoupled radial functions Z\ ' and Z\ ' are not simply degrees 

of freedom of the EM and gravitational field for any Q ^ 0. Only in the limit Q = do the potentials Vjf ^ and 
describe, respectively, purely EM and axial-gravitational perturbations of a Schwarzschild black hole. The radial 
equations (|f 3f) are solved by a series expansion of the form 



(-) 



(r+ - r_)~ 



L (r-r_) 



-/(r-4 



oo 



(16) 



where it = (r — r + )/(r — r_) and the coefficients Oj are determined by a four-term recursion relation. The problem 
can be reduced to a three- term recursion relation of the form (JSJ using a Gaussian elimination step j4ll | . Then we 
can use standard techniques to determine the quasinormal frequencies. The convergence of the summation can be 
improved using Nollert's expansion © for the rest. The first few coefficients of Nollert's series can be found in [30| . 



IV. FREEZING REISSNER-NORDSTROM: AN APPROXIMATE DECOUPLING 



If we are interested in the oscillations of an astrophysical black hole we can usually assume the black hole charge 
(in geometrized units) to be small, Q <§; M . In this limit it is reasonable to ignore the coupling between the EM field 
and the metric. In the following we study QNM frequencies of the RN black hole ignoring the mutual effect of the 
EM field on the metric and vice versa. We simply decouple the system Ijf 2(1 setting to zero the EM (gravitational) 

perturbations h[ ' (H2 )• The resulting equations are: 



dr 2 



V* 



H 



(-) 



0. 



(* = 1,2). 



(17) 



The gravitational potential when we freeze EM perturbations, h[ ^ = 0, is 

A 



|(| + l)r-3+^ 
r 



(18) 



V£ reduces to the standard RW potential for Q = 0. The EM potential freezing metric perturbations, H. 
becomes 



(-) _ 



r 



(19) 



Similarly to V^*, V* reduces to the potential for EM perturbations of Schwarzschild black holes when Q = 0. The two 
potentials can be written in the compact form 



I '* 



1(1 + l)r + (I - s 2 ) + 



4Q 2 



(i = |*| = 1,2). 



(20) 
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TABLE IV: Scalar, EM and gravitational QNMs of the RN black hole. The first line gives third order WKB results, and the 
second line sixth order WKB results for the decoupled potentials 12011 and 1211 . Sixth order WKB results have been kindly 
provided by Vitor Cardoso. The third line gives Leaver's results for the DF equation in the limit a — » 0. In the EM and 
gravitational cases, the fourth line gives results from Leaver's method applied to the coupled EM-gravitational potentials V/ - ' 
(i = 1,2). 



Scalar, 1 = 



= 



Q = 0.1 



= 0.2 



= 0.3 



Q = 0A 



WKB3, V 
WKB6, Vo 
Leaver, DF 



(0.209290 
(0.220939 
(0.220910 



-0.230390 
-0.201628 
■0.209791 



0.210984 
0.222496 
0.222479 



■0.230440) 
-0.202043) 
■0.210122) 



(0.216361 
(0.227572 
(0.227471 



■0.230392 
-0.203102 
0.210998 



0.226237 
0.237028 
0.236910 



-0.229299 
-0.204762 
■0.211849 



0.240493 
0.251727 
0.253168 



-0.223521) 
-0.207292) 
-0.210094) 



Scalar, I 



Q 



o 



Q 



o.i 



Q 



0.2 



Q 



0.3 



Q 



0.4 



WKB3, Vo 
WKB6, Vo 
Leaver, DF 



(0.966422 
(0.967284 
(0.967288 



-0.193610 
0.193532 
-0.193518 



0.973006 
0.973859 
0.973863 



-0.194022) 
-0.193948) 
-0.193933) 



(0.993994 
(0.994823 
(0.994826 



-0.195186 
-0.195120 
-0.195107 



1.033982 
1.034772 
1.034775 



-0.196707 
-0.196664 
-0.196649 



1.105332 
1.106104 
1.106105 



-0.196689) 
-0.196689) 
-0.196675) 



EM, I = 1 



= 



= 0.1 



= 0.2 



= 0.3 



Q = 0.4 



WKB3, Vr 
WKB6, VT 
Leaver, DF 
Leaver, V± 



(-) 



(0.491740 
(0.496383 
(0.496527 
(0.496527 



■0.186212 
■0.185274 
■0.184975 
-0.184975 



0.498236 
0.502809 
0.500367 
0.502950 



■0.186976) 
-0.186118) 
■0.185468) 
-0.185805) 



(0.519336 
(0.523726 
(0.512675 
(0.523843 



■0.189280 
-0.188648 
-0.186892 
-0.188312 



0.561324 
0.565304 
0.536404 
0.565513 



■0.193009 
-0.192902 
■0.188917 
-0.192408 



0.643293 
0.646616 
0.579729 
0.646987 



-0.196606) 
-0.197253) 
-0.189708) 
-0.196545) 



EM, I 



Q 







Q 



o.i 



Q 



0.2 



Q 



0.3 



Q = 0.4 



WKB3, VT 
WKB6, Vi* 
Leaver, DF 
Leaver, V± 



(-) 



(0.914262 
(0.915187 
(0.915191 
(0.915191 



-0.190130 
-0.190022 
0.190009 
-0.190009 



0.922188 
0.923104 
0.921657 
0.925930 



-0.190667) 
-0.190565) 
■0.190451) 
-0.190747) 



(0.947700 
(0.948588 
(0.942307 
(0.959852 



-0.192226 
-0.192142 
0.191710 
-0.192884 



0.997397 
0.998235 
0.981812 
1.024022 



-0.194478 
■0.194426 
0.193421 
-0.196033 



1.090249 
1.091011 
1.052923 
1.140260 



-0.195485) 
-0.195485) 
-0.193730) 
-0.198138) 



Gravitational, 1 = 2 



= 



= 0.1 



Q = 0.2 



Q = 0.3 



= 0.4 



WKB3, V 2 * 
WKB6, V 2 * 
Leaver, DF 



Leaver, V 5 



(-) 



(0.746324 
(0.747239 
(0.747343 
(0.747343 



-0.178435 
-0.177781 
0.177925 
-0.177925 



0.751578 
0.752520 
0.753643 
0.749489 



-0.178831) 
-0.178204) 
-0.178602) 
-0.178150) 



(0.768355 
(0.769379 
(0.773902 
(0.756874 



-0.179955 
-0.179415 
0.180620 
-0.178796 



0.800439 
0.801591 
0.813278 
0.772435 



-0.181460 
-0.181106 
■0.183793 
-0.179627 



0.858116 
0.859434 
0.886468 
0.802434 



-0.181608) 
-0.181619) 
-0.186603) 
-0.179286) 



For completeness we also consider scalar perturbations of the RN black hole: in this case no approximations are 
involved, and the potential reads 



A 

V = - s 



l(l + l)r + l 



2Q 



r 



(21) 



In Table ITVI we present the fundamental scalar, EM and gravitational QNM frequency for five selected values of the 
charge. From top to bottom we consider scalar perturbations with I = and I = 2, EM perturbations with I = 1 and 
I = 2 and gravitational perturbations with 1 = 2. 

We computed QNMs of the potentials l|20() and (|21|l using a third order WKB expansion then we checked 
convergence using a sixth order WKB expansion |45|. Sixth order results have been kindly provided by Vitor Cardoso. 
WKB results from the third (sixth) order expansion are given in the rows marked WKB3 (WKB6). "Leaver, DF" 
means that QNM frequencies have been computed using continued fractions and the a = limit of the DF equation 
PJl with the appropriate value of s. Finally, "Leaver, v/ means we applied Leaver's continued fraction method 
to the coupled EM-gravitational system Q13fl: the results are the "true" oscillation frequencies of the RN black hole 
I30L l4ll . Some comments are in order: 



i) The sixth order WKB results can be considered reliable. Even in the case in which we expect the worst 
convergence (scalar perturbations with I = 0) a sixth order WKB expansion agrees very well with results from 
the DF equation (which in this case, we stress it again, is exact). For scalar modes with I = 2 the agreement 
between the sixth order WKB and the results from the DF equation is quite astonishing. This is really a 
double check: not only it proves that WKB results can be considered reliable, it also shows (by a completely 
independent calculation) that the scalar QNMs obtained in Sec. Ull have the correct limit as a — > 0. 
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ii) The perfect agreement between decoupled and coupled EM perturbations with I — 1 is no surprise. Mathemat- 
ically, a quick inspection of Eq. (| 1 31) with i = 1 reveals that, since for I = 1 we have q-i = for any value of 
Q, EM perturbations with I — 1 are always decoupled from gravitational perturbations. Of course, physically 
this decoupling is due to the nonradiative character of dipolar gravitational fields (the first radiative multipole 
being the quadrupole, 1 = 2). Notice however that for EM fields and I — 1 QNMs of the DF equation disagree 
with the other two approaches, the deviations increasing for large values of Q. This is evidence that when we 
consider the DF equation with s = — 1 we are not really "killing" gravitational perturbations. However, to a good 
approximation, QNM frequencies of the DF equation for s = — 1 are very close to frequencies obtained freezing 
metric perturbations in RN when we consider small values of the charge, say Q < M/2. This is also true for 
EM QNM frequencies with I = 2, in which case EM perturbations are coupled with gravitational perturbations. 

iii) Inspection of gravitational QNMs with I = 2 confirms the above conclusions: when we consider the DF equation 
with s = —2 we are not really "killing" EM perturbations. However, to a good approximation, QNM frequencies 
of the DF equation for s = —2 are very close to those obtained freezing EM perturbations in RN when we 
consider small values of the charge, say Q < M/2. 

iv) For RN black holes with Q < M/2 the DF equation provides the correct RN oscillation frequencies of the full 
coupled EM-gravitational perturbations system within about 1 %. We can reasonably expect our results to have 
the same level of accuracy for rotating, KN black holes, at least when Q < M/2 and a is not too large. 



A. Asymptotic modes of "frozen" RN black holes 



Highly damped black hole QNMs received considerable attention recently, due to a conjectured relation with 
quantum gravity 0] . Motl and Neitzke ^3 used a monodromy calculation to compute asymptotic QNM frequencies 
of RN black holes. The key element of their calculation is the leading asymptotic behavior of the potential close to 
the origin. They showed that, since the leading term of both potentials in (|12|) as r — > is 

vs~4ir' (i=1,2) ' (22) 

with j — 5/3, asymptotic QNMs are given by the implicit formula 

eJ 3uJ + 2 + 3e- 0IU = . (23) 

where (3 and /?/ are the Hawking temperatures of the outer and inner RN horizons, respectively |47| . The "frozen" 
potentials i|20[l have the same leading-order behavior: 

AO 4 4 



where we used the fact that r* ~ r 3 /(3Q 2 ) as r — > 0. Therefore, if we kill EM (or gravitational) perturbations 
keeping the charge in the background, the asymptotic QNMs are still given by (|23|l . In this sense, the wild oscillations 
of asymptotic QNM frequencies as a function of Q are not induced by the EM-gravitational coupling. In fact, the 
presence of the charge term in the potential (and not the coupling of EM-gravitational perturbations) is the reason 
for the different topology of the Stokes lines in the monodromy calculation [48|. 



V. CONCLUSIONS 



In this paper we presented the first continued-fraction calculation of QNMs of the KN black hole based on the 
approximate perturbation equation derived by Dudley and Finley (16j. For scalar perturbations their equation is 
exact; for EM and gravitational perturbations, it provides a good approximation for values of the charge Q < M/2. 
We found no evidence for instabilities, extending a previous analysis by Mashhoon [20|. who presented arguments 
in favor of the stability of KN black holes in the eikonal approximation. To understand the meaning of the DF 
approximation we analysed in detail the zero-rotation (RN) limit. We computed QNMs for "frozen" (purely EM or 
gravitational) perturbations of RN black holes and compared results with the zero-rotation limit of the DF equation. 
We found that the DF equation is equivalent to frozen EM (gravitational) perturbations of RN only when Q < M/2. 
In this regime, and for a = 0, the DF equation provides the correct QNMs of the full coupled EM-gravitational 
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perturbations system within about 1 %. We expect our results to have the same level of accuracy also for rotating 
black holes, at least when Q < M/2 and a is not too large. 

To confirm this expectation we need a more complete analysis of the KN perturbation problem. Some su gge stions to 
decouple EM and gravitational perturbations of KN black holes were given in and, more recently, in |50j |. Chitre 
|5l| found a separable equation for rotating black holes with small values of Q. His treatment was later extended by 
Lee ( see a l so A s far as we know, no attempt has been made to compute the corresponding QNMs. Results 

obtained from Lee's small-Q expansion should not deviate much from ours, since the DF equation gives accurate 
QNM frequencies in the limit Q <C M. 

Our study shows that the DF equation can be used to provide reliable estimates of the QNM frequencies of slightly 
charged, rotating BHs. Mashhoon's approximate treatment, summarized in the Appendix, captures the essential 
physics and gives very good estimates for the I = \m\ = 2 fundamental QNM as a function of both Q and a (see 
Fig. |3) . These results could be useful to investigate gravitational wave emission in coincidence with black hole models 
of gamma-ray bursts. Proposed models to date include bipolar outflow from a fast-rotating, slightly charged KN black 
hole surrounded by a magnetosphere Q, black hole- torus systems P a i r creation due to vacuum polarization in 
gravitational collapse to a charged black hole ^lj and hyper-accretion scenarios 

An interesting extension of our work concerns the study of Dirac perturbations of the KN metric. In this case 
separability is not an issue |l5j |. Pekeris and collaborators considered the nucleus as a Kerr-Newman source 0. 
Assuming the angular momentum of the source to be the intrinsic spin angular momentum of the nucleus they found 
a remarkable result: the Dirac equation in the KN background predicts the hyperfine splitting observed in muonium, 
positronium and hydrogen to within the uncertainty in the respective QED corrections, except for an unaccounted 
factor of 2. Their analysis is worth revisiting, considering also that a very complete mathematical study of the Dirac 
equation in the KN spacetime appeared after Pekeris' work [54[. We hope to return to this problem in the future. 
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APPENDIX A: PERTURBED ORBITS OF NULL RAYS AND THE EIKONAL LIMIT 



In this Appendix we summarize Mashhoon's analysis of the perturbations of unstable circular photon orbits in the 
KN spacetime [l^. His physical model is closely related to (but slightly different from) Goebel's j2l|. The key idea 
is that, since QNMs represent a general property of the spacetime, we can find them using the time evolution of 
any convenient perturbation. In particular, Mashhoon considers an aggregate of massless particles in the unstable 
equatorial circular orbit of a KN black hole. The radius tq of this orbit is given by the roots of 



r 2 



- 3r /2 + 2Q 2 ± 2a(r /2 - Q 2 ) 1/2 = . (Al) 



The upper sign refers to corotating and the lower sign to counterrotating orbits, and only solutions with ro > r + are 
acceptable. The corresponding orbital frequencies are 

lJ± = a±rUrJ-Q^- (A2) 

According to Mashhoon's model, the real part of QNM frequencies in the eikonal limit is proportional to the frequency 
of the perturbed bundle of null rays that escape to infinity, w± , and the imaginary part is proportional to 

r (2r - 1) 

Perturbed null rays corotating (counterrotating) with the black hole correspond to I — ±m, so Mashhoon's analysis 
predicts QNM frequencies 

= (±liu ± ,-(n- l/2) 7 ), (A4) 

where the —1/2 in the imaginary part accounts for the fact that we count modes starting from n = 1, following 
Leaver's convention [2j|. In the Schwarzschild limit Mashhoon's model predicts ui± = 7 = 2/(3-\/3) ~ 0.384900, in 
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agreement with the leading-order large-/ expansion from the WKB approximation (Eq. (3.1) in [38 

(21 2(n- 1/2)" 



V3V3' 3V3 J ' ^ 

In the extremal limit we can eliminate (say) Q in favour of a in Eq. IjAlfl and solve analytically for the extremal 
frequencies Wj_ xt . For corotating orbits the result is w+ xt = (2 — 3a) _1 for < a < 1/4, and uff* = 4a(l + 4a 2 )" 1 for 
1/4 < a < 1/2. For counterrotating orbits, aj1 xt = —(2 + 3a) _1 in the whole range < a < 1/2. 
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